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CN Abstract 

O We analyze the influence of finite ion size effects in the response of a salt-free concentrated suspension 

^ of spherical particles to an oscillating electric field. Salt-free suspensions are just composed of charged 

£P colloidal particles and the added counterions released by the particles to the solution, that counterbal- 

^ ance their surface charge. In the frequency domain, we study the dynamic electrophoretic mobility of 

the particles and the dielectric response of the suspension. We find that the Maxwell- Wagner-O'Konski 

process associated with the counterions condensation layer, is enhanced for moderate to high particle 

charges, yielding an increment of the mobility for such frequencies. We also find that the increment of 

the mobility grows with ion size and particle charge. All these facts show the importance of including 

ion size effects in any extension attempting to improve standard electrokinetic models. 
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In the last years, there has been a renewed in- 
terest in electrophoresis. This is in part due to 
recent advances in nanoscience which make pos- 
sible the separation of macromolecules by size or 
charge. Suspended DNA or proteins are driven 
and separated by applying dc or ac electric fields 
[IH3]. Another main application in the field of 
nanoscience is the use of gold nanoparticles for 
drug delivery or cancer cell detection [U E] • Mea- 
surements of the electrophoretic mobility have 
been found to be useful to characterize the sur- 
face functionalization of these gold nanoparticles 
[6]. Usually particles are charged and suspended 
together with microions and a structure of elec- 
tric double layer (EDL) appears around the par- 
ticle 0IE]- The electrophoretic mobility of a sus- 
pended particle is not only dependent on the par- 
ticle charge or the viscosity of the medium, but 
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also on the configuration of the EDL. 

Most of the theoretical EDL models are based 
on the classical Poisson-Boltzmann equation, a 
mean field approach that takes into account point- 
like ions in solution. This theory breaks down 
when the crowding of ions becomes significant, 
and steric repulsion and correlations potentially 
become important. Some authors have shown 
that the consideration of finite ion size effects al- 
lows for the crowding of ions near the particle 
surface [9]. This redistribution of ions modifies 
the EDL around the particle and consequently its 
electrophoretic mobility when an external electric 
field is applied [TUtCE]. 

We can find in the literature different studies 
dealing with ion size effects. Some of them con- 
cern microscopic descriptions of ion-ion correla- 
tions [IH [T3]. These approaches are mainly re- 
stricted to equilibrium conditions, but are able 
to predict important phenomena like overcharg- 
ing [16J. Other studies are based on macroscopic 
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descriptions considering average interactions by 
mean field approximations [101 ITTH20] . In many 
of these works, the finite ion size is commonly 
included by modifying the activity coefficient of 
the ions in the electrochemical potential or by in- 
corporating entropic contributions related to the 
excluded volume of the ions. The macroscopic 
approaches have been found to work appreciably 
well with monovalent electrolytes for high particle 
charges and/or large ionic sizes when they have 
been compared with some simulation results [2T] . 
Most works in electrokinetics concern suspen- 
sions with low particle concentration, but nowa- 
days it is the concentrated regime that deserves 
more attention because of its practical applica- 
tions. These systems are difficult to understand 
due to the inherent complexity associated with 
the increasing particle-particle electrohydrody- 
namic interactions as particle concentration grows 
[221 123] ■ O n the other hand, systems with low salt 
concentration show a lower screening of the re- 
pulsive electrostatic particle-particle interactions, 
which favors the generation of colloidal crystals 
or glasses. Suspensions just composed of charged 
particles and their ionic countercharges (the so- 
called added counterions) in the liquid medium 
are named salt-free suspensions. The interest in 
these systems has increased in recent years from 
both experimental and theoretical points of view 



The dielectric response of colloidal suspensions 
as a function of the frequency of the applied ex- 
ternal electric field is a powerful tool. Its anal- 
ysis provides rich information on the dynamics 
of the EDL because it is very sensitive to the 
particle-solution interface. Carrique et al. have 
studied the dynamic properties of salt-free con- 
centrated suspensions with point-like ions when 
an ac electric field is applied [23 ED]- To our 
knowledge, the only theoretical work in the fre- 
quency domain considering ion size effects is the 
one of Aranda-Rascon et al. for dilute suspensions 
with electrolytes [32]. The same authors have 
shown that the consideration of a minimum ap- 
proach distance of ions to the particle surface, not 
necessarily equal to their effective radius in the 
bulk solution, can predict overcharging for high 



electrolyte concentrations and counterion valence 

Our aim in this paper is to analyze the influ- 
ence of finite ion size effects in the response of a 
salt-free concentrated suspension to an oscillating 
electric field. We will study specially the dynamic 
electrophoretic mobility of the particles and the 
dielectric response of the suspension in the fre- 
quency domain. Recently, we have studied the 
equilibrium EDL [TjJ] and the response to a static 
electric field [13] of this kind of suspensions with 
ion size effects. In this paper, we will extend our 
previous works to ac electric fields following the 
treatment of Carrique et al. for salt-free concen- 
trated suspensions with point-like ions [2H1 ED] ■ 

The plan of the paper is as follows. In Sec- 
tion [2] we describe the electrokinetic model to ac- 
count for ion size effects. We give details of the 
resolution method of the electrokinetic equations 
and define the quantities we calculate in Section 
[U The results of the numerical calculations are 
shown in Section [4] and analyzed upon changing 
particle surface charge density, particle volume 
fraction, and size of the counterions. In order to 
show the realm of the finite ion size effect in salt- 
free suspensions, the results are compared with 
standard predictions for point-like ions. Conclu- 
sions are presented in Section [5j 

2. Model 

2.1. Electrokinetic equations 

We use a cell model [33, 34J to study the macro- 
scopic properties of the suspension from appropri- 
ate averages of local properties in a representative 
cell. In this approach, each spherical particle of 
radius a is surrounded by a concentric shell of the 
liquid medium, having an outer radius h such that 
the particle/cell volume ratio in the cell is equal to 
the particle volume fraction throughout the entire 
suspension, that is 



(1) 



In this approximation, we simulate the hydro- 
dynamic and electrical interactions between par- 
ticles in the suspension by proper specification of 



boundary conditions at the outer surface of the 
cell. 

Let us consider a spherical charged particle of 
radius a, surface charge density a, mass density p p 
and relative permittivity e rp immersed in a salt- 
free medium of relative permittivity e rs , mass den- 
sity p s and viscosity rj, with only the presence of 
the added counterions of valence z c and drag co- 
efficient A c . We consider finite size counterions as 
spheres of radius R with a point charge at their 
center. By applying to the system an oscillat- 
ing electric field E e~ tu,t of angular frequency u, 
the particle moves with a velocity v e e~ luJt , the dy- 
namic electrophoretic velocity. The axes of the 
spherical coordinate system (r, 9, p) are fixed 
at the center of the particle, with the polar axis 
(9 = 0) parallel to the electric field. The solution 
of the problem at time t requires the knowledge, 
at every point r of the system, of the electric po- 
tential, \l/(r, t), the number density of counterions, 
n c (r,t), their drift velocity, v c (r,£), the fluid ve- 
locity, v(r, £), and the pressure, P(r, t). The elec- 
trokinetic equations connecting them are [7J 135] : 
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As we are interested in studying the linear re- 
sponse of the system to an electric field, we apply 
a perturbation scheme. Thus, each quantity X is 
written as the sum of its equilibrium value, X°, 
plus a perturbation term, SX, linearly dependent 
with the field multiplied by the term e~ lult , that 
represents the time dependent sinusoidal response 
of the stationary state [29J . 

We introduce the finite size of the counterions 
by considering their excluded volume and includ- 
ing the entropy of the solvent molecules in the free 
energy of the suspension, F = U — TS [H] 
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being n™ ax the maximum possible concentration 
of counterions due to the excluded volume effect, 
defined as n™ ax = V' 1 , where V is the average 
volume occupied by an ion in the solution. The 
last term in Eq. (J8J) is the one that accounts for 
the ion size effect, and was proposed earlier by 
Borukhov et al. [TS]. Performing the variations 



n c (r, £)v c (r, t) = n c {r, i)v(r, t) 



n c (r, t)] (4) of the free energy with respect to \E f0 (r) and n°(r), 
combining both resulting expressions and apply- 
ing spherical symmetry, we obtain 
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In these equations, p c (r,t) is the electrochem- 
ical potential of the counterions, eo is the vac- 
uum permittivity and e is the elementary electric 
charge. The drag coefficient A c is related to the 
diffusion coefficient by A c = ksT/D c , where ks 
is Boltzmann's constant, and T is the absolute 
temperature. 



tO^rs i _|_ 



cxp 






1 



(9) 



where b c is an unknown coefficient that represents 
the ionic concentration where the equilibrium 
electric potential is chosen to be zero. Details 
about this modified Poisson-Boltzmann equation 
including ion size effects can be found in Ref . [19] . 



To obtain the perturbation terms of the quanti- 
ties of interest, due to the symmetry of the prob- 
lem, we make use of the following spherical func- 
tions: h(r), (p c ( r ), and Y{r) [3B] 

v(r) = (v r ,i)j,^) = 

-h(r)Ecos6, - — (rh(r))Esm6,0 ) (10) 
r r dr J 

8fi c {r) = — z c e<p c (r)E cos 9 (11) 

<Jtf(r) = -Y(r)Ecos9 (12) 

with E = |E|. 

Substituting the above mentioned perturbation 
scheme into the differential electrokinetic equa- 
tions, Eqs. ([2])- ([6]), neglecting nonlinear pertur- 
bations terms, and making use of the symmetry 
conditions of the problem we obtain 
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where the C operator is defined by 



(14) 



(r) - Y{r)) (15) 



c 



d 2 2d 
dr 2 r dr 



(16) 



In the case of a static electric field, u = 0, Eqs. 



( 13 )-( 15 ) turn into the expressions obtained in 
Ref. [13]. For point-like counterions, n™ aa; = 00, 
these equations become those of Refs. J29J [30] . 

According to Ref. [10J, we incorporate a dis- 
tance of closest approach of the counterions to 
the particle surface, resulting from their finite 
size. We assume that the counterions cannot 
come closer to the surface of the particle than 
their effective hydration radius, R, and, therefore, 
the ionic concentration will be zero in the region 



between the particle surface, r = a, and the spher- 
ical surface, r = a + R, defined by the counterion 
effective radius. This reasoning implies that coun- 
terions are considered as spheres of radius R with 
a point charge at their center. 

With this consideration, we solve the electroki- 
netic equations, Eqs. (J9|, ( 13 )-( 15 ), only between 
r = a + R and r = b. When we address the prob- 
lem in the region between r = a and r = a+R, the 
equations to solve turn into the Laplace equation 
for the equilibrium electric potential, and equa- 
tions C(Ch(r)) = 0, <j> e (r) = 0, and CY(r) = for 
the rest of the spherical functions, because this re- 
gion is free of charge. We call FIS+L model this 
complete model that includes ion size effects and 
also considers the distance of closest approach of 
the counterions to the charged particle surface. 

2.2. Boundary conditions 

The boundary conditions needed to solve the 
electrokinetic equations are analogous, but deal- 
ing with complex quantities, to those described in 
Ref. [13] (Section II B) for the response of a salt- 
free concentrated suspension to a static electric 
field including ion size effects. At the particle sur- 
face, we apply the continuity of the electric poten- 
tial, the discontinuity of the normal component 
of the displacement vector, the non-slip condition 
for the fluid and the impenetrability of ions to 
the solid surface. On the outer surface of the cell, 
we use the Kuwabara's boundary conditions for 
the fluid velocity field and the Shilov-Zharkikh- 
Borkovskaya conditions for the perturbed electric 
potential. Finally, if we consider a distance of 
closest approach of the counterions to the particle 
surface, we also need the continuity of the pres- 
sure and both normal and tangential components 
of the fluid velocity and the vorticity. A remark- 
able difference with the static case is that the net 
force acting on the particle or the unit cell is not 
zero. Details of the net force calculation can be 
found in Ref. [32] or Appendix 1 in Ref. [29J . 



In terms of the radial functions \& (r), Y(r) 
4> c {r) and h(r), the boundary conditions are: 
(i) at the particle surface r = a 
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(18) 
(19) 
(20) 



(ii) at the surface r = a + R defined by the coun- 
terion effective radius 
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(iii) and finally, at the outer surface of the cell 



*°(6) = 
d^°(r) 



r=b 



ri— \rCh(i 
dr L 



r=b 



dr 
Y(b) = b 

Mb) = b 
Ch{b) = 

- z c eb c Y(b) = iup s 



(30) 

(31) 

(32) 
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This last boundary condition, Eq. (35), stands 



for the equation of motion of the unit cell. In the 
case of a static electric field, uj = 0, we recover the 
expression for the net force showed in Ref. |13j . 



3. Method and calculated quantities 

3.1. Method 

We will discuss the results of the proposed 
FIS+L electrokinetic model. In order to show the 
realm of the finite ion size effect in salt-free sus- 
pensions, the results are compared with standard 
predictions for point-like ions, PL model [29| 130] . 
The electrokinetic equations with their boundary 
conditions form a boundary value problem that 
can be solved numerically using the MATLAB 
routine bvp4c [37] . 

For the sake of simplicity, we assume that the 
average volume occupied by a counterion is V = 
(2R) 3 , being 1R the counterion effective diame- 
ter. With this consideration, the maximum pos- 
sible concentration of counterions due to the ex- 
cluded volume effect is n™ ax = (2R)~ 3 . This cor- 
responds to a simple cubic package (52% pack- 
ing). In molar concentrations, the values used in 
the calculations, n™ ax = 22, 4 and 1.7 M, cor- 
respond approximately to counterion effective di- 
ameters of 1R = 0.425, 0.75 and 1 nm, respec- 
tively. These are typical hydrated ionic diameters 
[38J . We present in Table [I] the parameter val- 
ues used in all the calculations. The chosen pa- 
rameters correspond to hydrated H + counterions, 
which are commonly found in many experimental 
conditions with salt-free suspensions of, for exam- 
ple, negatively charged sulfonated polymer parti- 
cles, due to the cleaning process of the suspension 
with proton exchange resins. 

3.2. Calculated quantities 

The dynamic electrophoretic mobility \x of a 
spherical particle in a concentrated colloidal sus- 
pension can be defined from the relation between 
the electrophoretic velocity of the particle and 



Table 1: Parameter values used in the calculations. 
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the macroscopic electric field. According to Refs. 
[T3l 129] it can be determined through 
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We calculate the nondimensional dynamic elec- 
trophoretic mobility as 
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The complex conductivity, K, of the suspension 
is usually defined in terms of the volume averages 
of the local electric current density and electric 
field in a cell representing the whole suspension. 
Following a similar procedure to that described 
for the dc conductivity in Ref. [13J, we obtain 
(see also Ref. [30] ) 
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From the complex conductivity, the real e' r (uj) 
and imaginary e"(ui) components of the complex 
relative permittivity of the suspension e r (u) are 
calculated by writing 

K{uj) = K(u> = 0) — iueoe r (co) 

= K(u = 0) + tue e"(u) - icoe e' r (cu) (39) 
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4. Results and discussion 

4-1. Point- like model 

The classical frequency response of a salt-free 
concentrated suspension with point-like counteri- 
ons, PL model, is as follows: (i) at low frequency, 
the electromigration and diffusion processes have 
enough time to be fully developed around the par- 
ticle and, commonly, this fact leads to the gener- 
ation of an induced electric dipole moment that 



tends to brake the particle motion. In this fre- 
quency region, there is a plateau value of the 
dynamic electrophoretic mobility that coincides 
with the electrophoretic mobility in static elec- 
tric fields, (ii) As the frequency increases, we 
find a frequency region where the counterions can- 
not follow the comparatively fast field oscillations. 
Thus, the above mentioned dipolar moment de- 
creases and, consequently, the dynamic mobility 
increases. This process is known as Maxwell- 
Wagner- O'Konski (MWO) relaxation and takes 
place whenever the medium and the charged par- 
ticle, surrounded by its EDL, present different 
conductivities and permittivities, (iii) Finally, the 
frequency can be so high that the inertia of the 
particle and fluid restricts the motion progres- 
sively. As a result, the mobility shows a contin- 
uous decline when the frequency rises, which is 
known as the inertial relaxation. 

There is another classical relaxation mech- 
anism, the alpha relaxation [7], which is re- 
lated to the concentration polarization effect (i.e., 
the presence of a gradient of neutral electrolyte 
around the particle). We do not find any alpha 
relaxation in a salt-free suspension, as was ex- 
plained in Ref. jlQ], because we only have one 
ionic species, the added counterions. 

In order to clarify the discussion, we will try to 
separate the different mechanisms by performing 
three different variations of the PL model: 

(i) the complete PL model that includes all the 
mentioned effects; 

(ii) the pure inertial response where we do not 
allow any perturbation of the ionic atmosphere 
from the equilibrium values. So, breaking mecha- 
nisms associated with the charge polarization are 
excluded in this variation; 

(iii) the inertia-free response, where we have 
eliminated all the inertial terms in the electroki- 
netic equations. 

In Fig. [I] we show the modulus of the scaled 
dynamic mobility for the three last-mentioned PL 
variations. We use two different particle volume 
fractions at a given particle charge density. We 
display in solid lines the complete PL model, in 
dashed lines the pure inertial response, and in dot- 
ted lines the inertia-free response. 
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Figure 1: Modulus of the scaled dynamic electrophoretic 
mobility as a function of frequency for different particle 
volume fractions. We use three variations of the PL model. 
Solid lines stand for results of the complete PL model. 
Dashed lines show the pure inertial response. Dotted lines 
display the inertia-free response. 



The pure inertial response behaves as follows: 
after an initial low frequency mobility plateau, the 
mobility monotonously decreases with frequency. 
This plateau has larger values than that of the 
complete PL model. The difference between the 
pure inertial response and the PL model is due 
to the absence of breaking effects on particle mo- 
tion associated with the induced dipole moment 
(double layer relaxation effect). 

In the numerical results corresponding to the 
inertia-free response, we observe one or two suc- 
cessive increments in the dynamic mobility for 
high or low particle volume fraction, respectively. 
These increments are related to one or two suc- 
cessive MWO relaxation processes in each case. 
When we introduce the inertia (dashed lines) to 
get the complete PL model (solid lines), these 
above mentioned increments become into one or 
two successive peaks, as can be observed. The 
diminution of the dynamic mobility after the first 
maximum is therefore due to the inertial relax- 
ation. 

As was expected, once all the MWO processes 
have relaxed in the high frequency limit, the mo- 
bility in the inertia-free response curve reaches 
the same plateau value than the one of the pure 
inertial response for low frequencies. This is be- 



cause, when the induced polarization completely 
disappears, the counterions distribution coincides 
with that of the equilibrium as in the pure inertial 
response. 

4-2. Maxwell- Wagner- O'Konski relaxations 

We have observed that two differenced MWO 
relaxations can exist. They will be related to two 
differenced regions in the EDL. Through the Wag- 
ner formula for a constant dielectric mixture, it is 
possible to obtain the frequency and the dielectric 
increment of a MWO relaxation [7] 
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where K p and K s are the conductivities of the par- 
ticles and the medium, respectively. The particles 
are assumed to be made of a nonconducting mate- 
rial. Their conductivity is exclusively associated 
with the surface conductivity, K a , that appears 
due to an excess of counterions in the ionic atmo- 
sphere, K p = 2K a /a 0. Eqs. (42b and (pb were 



derived without allowance of mutual polarization 
of particles and they are valid for suspensions with 
added electrolyte, thin EDL, and reasonably low 

The latter equations predict only one MWO 
relaxation process. However, for suspensions of 
highly charged particles, two different MWO re- 
laxations have been considered in the literature to 
explain their dielectric response. This considera- 
tion is based on the existence of two differenced 
regions in the EDL, specially when a is sufficiently 
high and ion size effects are considered, see Fig. 
2 on Ref. [19]. The first one is a condensate of 
counterions very close to the particle surface. The 
second one is a diffuse layer that extends from 
the end of the condensate to the outer surface of 
the cell. When we have finite size counterions, 
the condensate consists in a homogeneous region 



where counterions are well packaged. In the PL 
case, such picture of the condensate is not valid, 
but there are theoretical evidences of the existence 
of a thin region with different behavior in the elec- 
tric potential and ionic distribution than in the 
diffuse layer (39J SO]. As was suggested on Ref. 
[40], we will consider each region with a different 
associated MWO relaxation process and roughly 
estimate their MWO relaxation frequencies and 
dielectric increments to qualitatively explain the 
behavior of both MWO relaxation processes. 

To obtain the relaxation frequency of the con- 
densate, we need to calculate the surface conduc- 
tivity of the counterion condensation layer. Con- 
sidering this layer with a mean concentration ri % 
and a thickness 5, we obtain 
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For the study of the condensate relaxation, 
the conductivity of the counterions in the diffuse 
layer, K s , will be taken equal to zero, because it 
has no influence in the condensate relaxation pro- 
cess. Using K p = 2K a /a and e rp <C e rs , according 



to Eq. (42) we obtain a relaxation frequency 



to 
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and with Eq. (43) a dielectric increment 
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To obtain only the relaxation process of the dif- 
fuse layer, we need to calculate the conductivity 
of the counterions in this region, K s , 



K„ 



z 2 e 2 n r 



Xr 
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(47) 



where n c is the average counterions concentration 
in the diffuse layer, and a dl * is the charge den- 
sity at the spherical surface, r = a + 6, defined 
by the thickness of the condensation layer. We 
now take K p = because we consider the particle 
with its condensate as an equivalent particle with 
less surface charge density, a * . Introducing the 



expression of the conductivity K s in Eq. (42) and 
considering e rp <C e rs , we obtain the relaxation 
frequency of the diffuse layer 



L) 



dif 
MWO 



— 3z c ea dt f (j) 
e e rs aX c (l - 4 



(48) 



and with Eq. (43) the dielectric increment 



A */ _ 90(1 - 0K 2 P 
MWO " (2 + 0)3e rs 



(49) 



As commented before, the previous expressions 
will be more precise for reasonably low 0. Also, 
the expression for the MWO relaxation frequency 
of the condensate will be more accurate than the 
one for the diffuse layer: the use of an average 
concentration works better in the condensate be- 
cause it is a thin layer with homogeneous ionic 
density. 

4-3. Finite ion size 

We will jointly study both the dynamic elec- 
trophoretic mobility of the particles and the di- 
electric response of the suspension as a function of 
frequency, because they are strongly interrelated. 

4-3.1. Condensate MWO relaxation 

Fig. [2] represents the modulus of the scaled dy- 
namic electrophoretic mobility and Fig. [3] the real 
(a) and the imaginary (b) parts of the relative per- 
mittivity of a salt-free concentrated suspension as 
a function of frequency. We compare the results of 
the FIS+L model (dashed lines) at a fixed counte- 
rion size, n™ ax = 4 M, with those of the PL model 
(solid lines). Different colors stand for different 
particle surface charges. The calculations were 
made at low particle volume fraction, = 10~ 2 . 

For the low frequency mobility and permittiv- 
ity plateaus, Figs. [2]and[3^, we observe that there 
is almost no difference between the finite size and 
point-like results for any particle charge at low 
particle volume fraction. As was discussed in Ref. 
[13j for the static mobility and conductivity, this 
is because in the case of a dilute suspension the 
inclusion of ion size effects only significantly mod- 
ifies the counterions fluxes in the immediate vicin- 
ity of the particle. 
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Figure 2: Modulus of the scaled dynamic electrophoretic 
mobility of the particles as a function of frequency for dif- 
ferent particle surface charge densities. All calculations 
performed at low particle volume fraction. Solid lines show 
the results for point-like ions. Dashed lines show the re- 



sults of the FIS+L model with 



= 4M. 



The MWO relaxation frequency is defined as 
that of the maximum in the imaginary part of the 
permittivity, as it is well-known. We note that for 
a suspension with low volume fraction, depending 
on particle surface charge, one or two differenced 
MWO relaxations, Fig. [3j or analogously one or 
two mobility maximums, Fig. [2j may take place. 
As we indicated before, we associate the first one 
upon increasing frequency with the MWO relax- 
ation of the diffuse part of the EDL, and the sec- 
ond one with the relaxation of the condensate. 

At low a there is no condensate of counterions 
near the particle surface and therefore no conden- 
sate MWO relaxation process is observed. When 
we rise the particle charge, almost all the extra 
counterions accumulate in the condensate [39|, 
which seriously grows when also ion size effects 
are considered [12] . This explains why ion size ef- 
fects do not produce any remarkable effect in the 
MWO relaxation of the diffuse part of the EDL 
and why they considerably enhance the conden- 
sate MWO relaxation. 

According to Eq. (48), the frequency of the 
MWO of the diffuse layer changes with <j dt f . As 
the counterions concentration in the diffuse layer 
has been scarcely altered, i^mwo remains basically 
the same in Figs. [2] and [3j On the contrary, Fig. 
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Figure 3: Real (a) and imaginary (b) parts of the relative 
permittivity of the suspension as a function of frequency 
for different particle surface charge densities. All calcula- 
tions performed at low particle volume fraction. Solid lines 
show the results for point-like ions. Dashed lines show the 
results of the FIS+L model with n™ ax = 4 M. 



^jp shows an increment of the frequency of the con- 
densate MWO relaxation when ion size effects are 



considered. This is in agreement with Eq. (45), 



because i^mwo increases when the width of the 
condensate, S, raises. For the well resolved MWO 
peaks of the two highest surface charge curves in 
Fig. [3]d, we observe that the height of the peaks 
is nearly independent on both, the surface charge 
density and the ion size, in accordance with Eqs. 



(46) and (49) 



In a previous paper, we studied the effects of 
the electric polarization on the magnitude of the 
static electrophoretic mobility in a salt-free con- 
centrated suspension with finite ion size effects 
[13J. This study was based on a procedure devel- 
oped by Bradshaw-Hajek et al. [42]. We showed 
that the induced charge polarization density was 
larger when ion size effects were considered. The 
generalization of the latter study to ac electric 
fields leads to similar conclusions. As the re- 
laxation effect and, correspondingly, the induced 
dipole moment, have been increased with ion size, 
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Figure 4: Modulus of the scaled dynamic electrophoretic 
mobility of the particles as a function of frequency for dif- 
ferent particle volume fractions. Solid lines show the re- 
sults for point-like ions. Dashed lines show the results of 
the FIS+L model with rC ax = 4 M. 
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Figure 5: Imaginary part of the relative permittivity of the 
suspension as a function of frequency for different particle 
volume fractions. Solid lines show the results for point-like 
ions. Dashed lines show the results of the FIS+L model 
with nT ax = 4 M. 



the second mobility maximum in Fig. [2] attains 
higher values as well. This is due to the disap- 
pearance of breaking mechanisms on the particle 
motion of increasing importance as the ion size 
grows, causing the dynamic mobility to reach su- 
perior values when the size of the counterions is 
taken into account. 

4-3.2. Overlapping of MWO relaxations 

We show the frequency response of the modu- 
lus of the scaled dynamic electrophoretic mobility 
and the imaginary part of the relative permittiv- 
ity of the suspension in Figs. [4] and |5j respec- 
tively. In both Figures we compare the results 
of the FIS+L model with n™ ax = 4 M, dashed 
lines, with those of the PL model, solid lines, for 
different particle volume fractions. All the calcu- 
lations were performed at a high particle charge, 
a = -40 /iC/cm 2 . 

We see in Figs. [2] and [3] how two differenced 
MWO relaxations take place when surface charge 
increases in conditions of low volume fraction. 
Now we observe that the MWO relaxations of 
the condensate and the diffuse layer in Fig. [5j or 
analogously the two mobility maximums in Fig. 
[4], tend to overlap in frequency for concentrated 
suspensions at high surface charge. According to 
Eq. (48), Co>M WO g rows with volume fraction at 



a rate 0/(1 — 0). This is in agreement with the 
frequency shift observed in the MWO relaxation 
of the diffuse layer, indicated with black arrows 
in Fig. [5j Eq. (45) predicts a frequency change 

0)/ (2 + 0) for 
Then we find 



with volume fraction at a rate (1 
the condensate MWO relaxation 



no significant changes in w^wo f° r l° w 4> values, 
and a small decrease for high volume fractions as 
shown in Fig. [5] These behaviors result in the 
observed overlapping of the MWO relaxations for 
concentrated suspensions. 

When we include ion size effects we find only 
changes in the condensate MWO relaxation (en- 
hancement of the corresponding mobility maxi- 
mum and small increase in w^wo)- These changes 
can be explained with the same reasoning used 
for Figs. [2] and [3] the consideration of finite 
size counterions seriously enlarges the condensate 
near the particle but does not produce remark- 
able effects in the diffuse layer. We also observe 
the well-known diminution of mobility with the 
increase of volume fraction in Fig. [4] due basi- 
cally to the larger screening of the particle charge 
[19J: when the particle concentration grows, the 
available space for the counterions inside the cell 
decreases and, consequently, the screening of the 
particle charge is greatly raised, thus reducing the 
value of the surface potential and, therefore, the 
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Figure 6: Modulus of the scaled dynamic electrophoretic 
mobility of the particles as a function of frequency for dif- 
ferent particle surface charge densities. All calculations 
performed at high particle volume fraction. Solid lines 
show the results for point-like ions. Dashed lines show the 
results of the FIS+L model with nT ax = 4 M. 



mobility. 

4-3.3. Highly charged concentrated suspensions 

Figs. [6] and [7] show the modulus of the scaled 
dynamic electrophoretic mobility and the real (a) 
and imaginary (b) parts of the relative permit- 
tivity of the suspension as a function of frequency 
for different particle surface charges. We compare 
the results of the FIS+L model with n™ ax = 4 M, 
dashed lines, with those of the PL model, solid 
lines at a high particle volume fraction, <fi = 0.5. 
Besides, Figs. [8] and [9] represent the same 
quantities at fixed particle surface charge, a = 
—40 /iC/cm 2 , and volume fraction, <fi = 0.5. 
In these Figures, we compare the results of the 
FIS+L model at different ion sizes (different col- 
ored dashed lines) with those of the PL model 
(black lines). 

As we mentioned before, for this high volume 
fraction value, the two MWO relaxations are over- 
lapped in a unique broad peak. We observe how 
both the dynamic mobility and relative permit- 
tivity increase when we consider finite size coun- 
terions in comparison with the PL case and when 
we increase the particle surface charge. The rea- 
son is that the increase of the surface charge or 
the consideration of ion size effects leads to an 
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Figure 7: Real (a) and imaginary (b) parts of the relative 
permittivity of the suspension as a function of frequency 
for different particle surface charge densities. All calcu- 
lations performed at high particle volume fraction. Solid 
lines show the results for point-like ions. Dashed lines show 
the results of the FIS+L model with n™ ax = 4 M. 



enhancement of the overall charge polarization in 
the EDL, resulting in both, higher permittivity 
values as can be seen in Figs. [7^, and [9^,, and 
larger heights of the corresponding peaks of the 
imaginary part, Figs. [TJd and[9]3. A similar expla- 
nation applies to the remarkable increment ob- 
served in the mobility maxima, associated with 
the MWO relaxations, Figs. [6] and |8j the dis- 
appearance of the augmented induced dipole mo- 
ment gives rise to greater mobility values. 

Figs. [7Jd and [9]b display a shift to larger fre- 
quencies in the MWO relaxation when ion size 
effects are considered. This shift is also larger 
the larger the size of the counterions (lower n™ ax 
value). The MWO relaxation of the diffuse layer is 



nearly independent of the ion size, Eq. (48), and, 



therefore, the shift observed is entirely due to the 
displacement of the condensate MWO relaxation. 
We have checked numerically that the expression 
of the MWO relaxation frequency of the conden- 



sate, Eq. (45), predicts a shift to higher frequen- 



cies with the increase of the ion size through the 
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Figure 8: Modulus of the scaled dynamic electrophoretic 
mobility of the particles as a function of frequency for dif- 
ferent ion sizes (dashed lines) . Black lines show the results 
for point like ions. 



product n™ ax 5. This is because the width of the 
condensate 5 augments in a higher rate than the 
parameter n™ ax diminishes upon increasing ion 
size [19] . When we increase the particle surface 
charge, we are also increasing the surface con- 
ductivity K u of the condensate, for both PL and 
FIS+L models, and consequently the frequency of 
the MWO relaxation rises. 



5. Conclusions 

By using a cell model approach we have an- 
alyzed the influence of finite ion size effects in 
the response of a salt-free concentrated suspen- 
sion of spherical particles to an oscillating electric 
field. We have derived a mean-field ac electroki- 
netic model that accounts for the excluded volume 
of the counterions. 

In the frequency domain, we have studied the 
dynamic electrophoretic mobility of the particles 
and the dielectric response of the suspension. For 
this purpose we have performed a comparative 
study of the different physical mechanisms, pure 
inertia response and charge polarization relax- 
ations, to know how they interplay to give the 
complete response. This study has allowed us 
to characterize the relative importance and relax- 
ation frequencies of each mechanism separately. 
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Figure 9: Real (a) and imaginary (b) parts of the relative 
permittivity of the suspension as a function of frequency 
for different ion sizes (dashed lines). Black lines show the 
results for point like ions. 



In the discussion of the numerical results two dif- 
ferent MWO relaxations have been successfully 
associated with the relaxations of the different 
ionic processes that take place in the diffuse and 
condensate regions of the EDL. Furthermore, the 
inclusion of ion size effects leads to an overall in- 
crement of the dynamic mobility and relative per- 
mittivity in comparison with the point-like case. 
The enhancement of the MWO relaxation for 
moderate to high particle charges, which is associ- 
ated with the counterions condensation layer, has 
yielded a remarkable increment of the mobility for 
such frequencies. In addition, we have found that 
this increment of the mobility grows with ion size 
and particle charge. Besides, we have observed a 
shift in the MWO relaxation of the condensate to 
larger frequencies with ion size. 

Some of these calculations can be compared 
with experimental results. To perform such 
comparisons, concentrated suspensions of highly 
charged particles are required. These suspen- 
sions have been classically difficult to synthe- 
size, although existing highly charged sulfonated 
polystyrene latexes are good candidates. More- 
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over, high-frequency experimental setups are 
needed to work with the very high frequency re- 
gion where the MWO relaxation of the condensate 
takes place. 
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